Quantum phases of bosons in double-well optical lattices 
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We study the superfluid to Mott insulator transition of bosons in a two-legged ladder optical lat- 
tice, of a type accessible in current experiments on double- well optical lattices. The zero- temperature 
phase diagram is mapped out, with a focus on its dependence upon interchain hopping and the tilt 
between double wells. We find that the unit-filling Mott phase exhibits a non-monotonic behavior 
as a function of the tilt parameter, producing a reentrant phase transition between Mott insulator 
and superfluid phases. 

PACS numbers: 03.75.Hh, 03.75.Lm, 05.30.Jp, 73.43.Nq 



Optical lattices loaded with ultra-cold atoms provide 
the opportunity to study quantum phases of many- 
particle systems because of their unprecedented degree 
of controllability [1]. Presently the lattice depth, dimen- 
sionality, geometry, and filling factor can all be reason- 
ably controlled. While one of the first examples of this 
degree of control was the experimental observation of 
the superfluid (SF)-to-Mott insulator (MI) transition in 
three-dimensional cubic optical lattices as a function of 
the lattice depth [2], tetragonal and orthorhombic opti- 
cal lattices can also be produced by deepening the optical 
potential along desired directions [3, 4]. 

More recently, possibilities for control have expanded 
with the experimental realization of double-well optical 
lattices. Control of the polarization of the laser beams 
allows for the production of lattices with a base in two 
and three dimensions [5] as illustrated in Fig. [1^, where 
Bose atoms (^^Rb) have been successfully trapped. In 
particular, one can create a one-dimensional double-well 
optical lattice corresponding to a two-leg ladder struc- 
ture shown in Fig. [TJd by increasing the long period of 
the double- well optical lattice. In standard condensed- 
matter physics, a few compounds, such as vanadyl 
pyrophosphate (VO)2P207 [y] and some cuprates like 
SrCu203 [7], have such two- leg ladders in their crys- 
talline structure, and they display much of the interesting 
physics encountered in general ladder systems, associ- 
ated with the interplay between spin-gapped and super- 
conducting states [8]. However, conventional condensed- 
matter systems come with fixed dynamical and structural 
parameters, while the flexible variability of optical lat- 
tices offers the prospect of exploring the full parameter 
space. Moreover, the particles confined in the current 
double-well optical lattices are bosonic atoms, in con- 
trast to conventional condensed-matter systems, where 
electrons (fermions) dictate the quantum phases. 

In this paper, we study the zero-temperature phase di- 
agram of bosons in double- well optical lattices. We focus 
on the case of a two-legged ladder, where analytical and 
numerical progress can be made; in particular, we ap- 




FIG. 1: (a) Schematic picture of a 3D configuration of a 
double- well lattice. Circles represent sites, and solid lines 
represent connections via the hopping between sites, (b) 
Schematic picture of the two-leg ladder, (c) Double-well po- 
tential corresponding to the cross section for the dashed line 
in (a) . Dotted lines represent the energy levels for each well. 



ply the time-evolving block decimation (TEBD) [l2|, [l^l 
method to such ladder systems. We show that the phase 
diagram changes dramatically as a function of the chem- 
ical potential /i, the intrachain (interchain) hopping t\\ 
(tx), the on-site repulsion U and the tilt A of the double- 
well, which are indicated in Figs. [TJd and [It. For A = 
and different ratios t^jU ^ Mott phases with half-integer 
(in addition to integer) filling factors emerge in the phase 
diagram oi fi/U versus tw/U for small ratios oit\\/U . As 
t^/U increases, the half- filling MI phase becomes larger 
and the unit-filling one becomes smaller. For fixed ra- 
tio t^/U ^ and different values of A//7, we also obtain the 
II /U versus t\\/U phase diagram which reveals a reentrant 
phase transition for the unit-filling MI induced by the tilt 
A. The reentrant phase transition can be attributed to 
the development of coherence in each double well in the 
vicinity oi \ = U which drives the system into the SF 
phase. Finally, we also calculate the critical points for 
the MI-to-SF transition at half and unit fillings. 

To describe the physics discussed above, we introduce 
the Bose-Hubbard model for the double-well ladder 



E 

ve{L,R} 



(1) 



2 



where Hf^ represents the double-weh Hamitonian for a 
given ladder index i and is given by 
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at creates a boson at the lowest level localized on the 
left (right) of the z-th double- well when r] = L (R). We 
do not include the effect of the harmonic trapping poten- 
tial. We assume that all the parameters are sufficiently 
small compared to the energy gap A between the first 
and second levels of each site. Furthermore, all param- 
eters of H are controllable in experiments [5[, and thus 
we begin our discussion by analyzing the limit of =0. 

Integer and Half-Integer Mott Phases: When = 
several MI phases emerge. In this case, the Hamilto- 
nian Eq. Q reduces to = Y^^Hf^. One can ob- 
tain the eigenenergy E^^\n^j) and eigenstate \^n,j) = 
'l2nL=o ^ni,(^7 j)|^L, n—riL) of Hf^ ^ where n is the num- 
ber of atoms in each double well and the quantum number 
j is a non- negative integer less than n + 1 (j = 0, 1, n). 
I^L, ^r) is the Fock state with ul (pr) atoms on the left 
(right) of each double well. When n is even (odd), the 
filling factor v of the MI phases is integer (half-integer). 
Although there exists a MI phase for each value of n, we 
focus here only on the v = 1, 1/2 phases. 

We consider first the case of symmetric double wells 
(A = 0) and discuss two limiting situations corresponding 
to ^ ^ and t± <C U. For t± /7, the anti-bonding 
single particle state of the double-well is pushed to en- 
ergies much higher than [/, and only the bonding single 
particle state and the lowest energy two-particle state are 
important. Therefore, in this limit, the u = 1/2 {u = 1) 
phase can be mapped into the unit (double)-filling MI 
phase for a ID lattice with an effective hopping ty and 
an on-site repulsive interaction U/2. For t_\_ <C /7, the 
u = 1 MI phase approaches that of a ID lattice (two in- 
dependent filling-one chains), and the width of the u = 1 
MI phase on the /i-line is ~ U. In the strict limit of 
= 0, the = 1/2 MI phase vanishes and the system is 
always a superfluid since there is a low energy path for 
the bosons to move along the chains. 

Next, we consider the case of tilted double wells, where 
A 7^ 0. The MI states present in the double- well ladders 
discussed here are quite different from those encountered 
in strictly ID superlattices Q. When A max(t^,/7), 
the MI with filling v in the double- well ladder reduces to 
the MI with 2iy in a single ID lattice. In this regime, a 
transition to a SF phase occurs at a crictical ty in contrast 
to the case of the ID superlattice, where all the occupied 
wells are completely isolated from each other and the 
system remains always in the MI phase. 

We discuss two special cases X = U with t± <^ U and 
A max(tx,t^). For X = U and t± <C t/, where the 
states |1, 1) and |0, 2) are nearly degenerate, the width of 
the 1^ = 1 MI phase on the /i-line is reduced to ~ 2V2t±. 



However, for A max(tx,^/), the width of the u = 1 
MI phase is ~ U. This happens because the u = 1/2, 1 
MI phases in a double-well ladder reduce to the unit- 
and double-filling MI phases of a single ID lattice, as all 
bosons prefer to be in the lowest energy side of the largely 
tilted double well. 

These special cases reflect the more general trend that 
as A increases, width of the u = 1 MI phase on the /i-line 
changes non-monotonically. Such non-monotonic behav- 
ior for the V = 1 MI phase is also found in (/i, ty )-plane 
for varying values of A and will be discussed next by tak- 
ing into account finite t\\ and studying the insulator to 
superfluid transition. 

Insulator to Superfluid Transition- 1: To include the ef- 
fects of ty and study the MI-to-SF transition, we use 
first a perturbative mean- field approach (PMFA) [T^. 
Althoiigh PMFA fails to describe ID systems quantita- 
tively [11] , it provides qualitative understanding and ana- 
lytical insight. Quantitative results can be obtained using 
the TEBD method 0, ^ ; these results are compared 
later with the picture that emerges from PMFA. 

We consider the effects of finite ty and introduce the 

SF order parameter ip^^ = {ai^r/) = {o^\r]) i^^^ the Hamil- 
tonian Eq. ([1]), which reduces to 

Hc,Y,Hf = J2[Ht- + 2t|| J2^l + n (3) 

i i 77 

where Vi = —2t\\ i^rjictl +<^i) describes the transfer of 
atoms between i-th sites and the condensate ?/^^ and is 
treated perturbatively. 

Using perturbation theory, we obtain the correction 
AEn = En — E^^"^ (n, 0) to the unperturbed ground state 
energy E^^\n^O) in terms of iprj- Performing a linear 
transformation (^1,^2)* = ^{iIjl^'^rY to diagonalize 
the quadratic part of AEn leads to 



AEn = Ac(n,f^,fy,/i)^^. 
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where X is a 2 by 2 Hermitian matrix and the bars on 
parameters mean the normalization by /7, e.g. Jl = ja/U. 
Expressions for the coefficients of and fourth order 
terms are quite long, thus we will not give them here. 
However, A2 is always positive, while Ai changes sign, 
and the fourth order coefficients are positive, leading to 
second-order phase transitions between the MI (^1 = 
^2 = 0) and SF (^1 7^ 0, ^2 = 0) states. 

For symmetric double wells (A = 0), we obtain analyt- 
ical expressions for the MI-SF phase boundary {Ai = 0) 
in two limits. When t± <C /7, the phase boundaries are 



, n = l (^ = 1/2), 
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In this case, the critical value of is 



n = l (^ = 1/2), 
U - \t^, n = 2 {v = 
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FIG. 2: Phase diagrams for symmetric double wells. Dashed 
lines represent the phase boundary calculated within PMFA. 
Solid lines, density plots, and dots are calculated by the 
infinite-TEBD. h is integer inside the solid lines. The den- 
sity plots represent a. The dots represent the critical point of 
the BKT transition. 



FIG. 4: Phase diagrams for tilted double wells (t± = O.lt/). 
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FIG. 3: Critical intra-chain hopping t|j for symmetric double 
wells as a function of t^. Dashed and solid lines represent 
tjj's for the half- and unit-filling MI phases calculated within 
PMFA. Triangles and squares represent t|j's for the half- and 
unit-filling MI phases calculated by the infinite-TEBD. 



with /ic — for = 1/2 and jic — (a/2 — for = 1. 
In the case of t±_ ^ [/, the double- well system reduces ef- 
fectively to a single ID lattice, and the phase boundaries 
as well as the value tjj can be obtained from the standard 
results [l^ by replacing U 

In Fig. [21 the MI-SF phase boundaries for A = calcu- 
lated within PMFA are shown as dashed lines for different 
values of t^. Notice that the figures are not in the same 
scale. In Fig. [3l the critical intra-chain hoppings tjj's 
for the V = 1/2, 1 MI phases are shown as functions of 
as dashed and solid lines. The v = 1/2 {v = 1) MI 
lobe grows (shrinks) since tjj increases (decreases) with 
increasing t±_ so that the double-well system reduces to 
a single ID system in the limit of ^ U . 

Next, we discuss the case of tilted double wells (A 7^ 0). 




X/U 

FIG. 5: t|j's for t± — 0.1 [/ as a function of A (Error bars are 
smaller than the size of symbols). 



In Fig.m we show the MI-SF phase boundaries for differ- 
ent values of A at fixed = 0.1 /7 indicated by dashed 
lines. In Fig.O we show the critical intra-chain hoppings 
tjj's versus A for the v = 1/2, 1 MI phases indicated by 
dashed and solid lines, respectively. The u = 1/2 MI 
lobe or tjj grows monotonically as A increases. In con- 
trast, the u = 1 MI lobe or tjj changes non- monotonically 
as a function of A. 

This non-monotonic behavior for = 1 can be under- 
stood as follows: at A = 0, tjj is given approximately by 
Eq. ([6]) since t± <^ U when t± =0.1 U. As A increases, 
tjj initially decreases. At A = /7, ty reaches a minimum 

when tjj ^/2t±/6, since the states |1, 1) and |0, 2) states 
are nearly degenerate, i.e., the local state of the MI phase 
at A = [/is |$2,o) ^ (|l,l) + |0,2))/v/2. The development 
of this local coherence then pushes the system into the SF 
phase. Further increase of A moves the system away from 
this degeneracy which favors SF, and forces tjj to increase 
causing a reentrance into a MI phase. In particular, when 
A [/ all atoms move to a single chain and are in |0, 2), 
thus the critical value becomes tjj :^ (5 — 2a/6)/7/2 as 



4 



expected for a single chain flO^. 

The non-monotonic behavior of tjj shows a reentrant 
quantum phase transition from MI to SF to MI, induced 
by the tilt A when tjj is kept between (tpmin and (tpmax- 
Taking into account the high degree of control achieved in 
double- well optical lattices [5] , we expect this reentrance 
to be observed experimentally. However, since we do not 
expect PMFA to give quantitatively correct results for 
the double-well (ladder) optical lattice, we next discuss 
numerical results using the TEBD method. 

Insulator to Super fluid Transition- II: To determine 
quantitatively the phase diagrams for double-well (lad- 
der) optical lattices, we use the infinite-size version of 
TEBD [13], which provides an excellent ground state 
for ID quantum lattice systems via imaginary time evo- 
lution. To apply the TEBD method to our problem, 
we map the double-well (ladder) Bose-Hubbard model 
into a single chain with nex^-neare^^-neighbor hopping, 
whose ground state can be calculated via the swap- 
ping technique We note that the infinite-TEBD 
algorithm has been recently applied to single chains 
with only nearest-neighbor hopping, where the quantum 
Berezinskii-Kosterlitz-Thouless (BKT) critical point [lH 
was obtained for the unit-filling case. While the maxi- 
mum number of bosons per site is iVmax = oo, conver- 
gence is already achieved in our numerical calculations, 
when A/'max = 4 for = 1/2 and A/max = 5 for = 1. 

The phase diagrams in the (/i, t||)-plane are shown in 
Figs. [2] and [H where the solid lines indicate the MI-SF 
phase boundaries, which have roughly a triangular shape. 
The sides of the MI lobe, the phase transition occurs from 
a 1/ = 1/2, 1 MI to a SF with 7^ 1/2, 1. However, the 
two sides of the "triangle" merge for each MI phase (see 
dots in Figs. 2 and 4) producing a phase transition from 
a z/ = 1/2, 1 MI to a SF with z^ = 1/2, 1, which is of the 
BKT type 

To locate the phase boundaries we calculate directly 
the mean number of atoms per double well h = ^^{ni^r])^ 

but we also calculate the fiuctuation a = (nf) — {fii)'^^ 
which is small deep in the MI regions, and relatively large 
in the SF regions. Since we are interested in local observ- 
ables, such as h and a, convergence is already achieved 
for X = 15, where x is the size of the basis set retained 
in the TEBD procedure [l2^. 

We locate the BKT transition on the lines of integer 
h {v = 1/2,1) by calculating the correlation function 
(a]^ao), where a\ creates an atom in the lowest single par- 



ticle state of a double-well. The SF phase of the double- 
well ladder can be regarded as a two-band Tomonaga- 
Luttinger liquid (TLL), and the correlation function ex- 
hibits power-law decay as (dj^ao) cx . The expo- 
nents i^c at the phase transitions can be calculated from 
the TLL theory. For instance, when max(t^,A) ^ [/, 
our system is effectively a single ID chain and has the 
critical value = 1/2 for the BKT transition [ll|. In 
addition, when A = and n = 2 {u = 1)^ the critical 
value is also Kc = 1/2 [16]. Consequently we use the 
criterion = 1/2 to identify the critical point for the 
BKT transition at integer values n (y = 1/2, 1). 

We calculate as a function of by fitting a ■ t~^I^ 
to the correlation function calculated from the TEBD 
method with x = 60. We use the intervals 10 < r < 15, 
15 < r < 20, 20 < r < 25, and 25 < r < 30 for the fit 
and take the average of them to produce error bars. The 
critical intra-chain hopping tjj along the lines of integer 
n is determined when K = Kc. The dots in Figs. [2] and 
m correspond to the BKT transition points. In Figs. [3] 
and [5l tjj's for n = 1,2 (i/ = 1/2,1) are shown as tri- 
angles and squares, respectively. The phase boundaries 
asymptotically approach those of PMFA as tends to 
zero. On the other hand, differences between PMFA and 
TEBD are significant when is relatively large. In par- 
ticular, t|j obtained using TEBD is more than twice as 
large as that obtained within PMFA. However, the qual- 
itative behavior of the phase diagram as a function of t± 
and A obtained within PMFA is consistent with that of 
the infinite-TEBD. 

Conclusions: In summary, we have studied the 
superfiuid-to-Mott insulator transition of bosons in 
double-well (ladder) optical lattices. Applying the time- 
evolving block decimation (TEBD) method to the two- 
leg Bose-Hubbard model, we have calculated the zero- 
temperature phase diagram. We have found that the 
phase diagram changes significantly depending on the 
inter-chain hopping and tilt of the double wells. In par- 
ticular, we have shown that the tilt can be used to in- 
duce reentrant transitions between Mott insulator and 
superfiuid phases. Through a comparison of the re- 
sults of TEBD and the perturbative mean-field approach 
(PMFA), we have shown that PMFA fails to describe the 
phase diagram quantitatively, but captures its qualitative 
trends. 
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